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o 

Abstract 

o ■ 

Outer envelopes of neutron stars consist mostly of fully ionized, strongly coupled Coulomb 
plasmas characterized by typical densities p ~ 10 4 -10 n gem -3 and temperatures T ~ 10 4 - 
10 9 K. Many neutron stars possess magnetic fields B ~ 10 11 -10 14 G. Recent theoretical 
advances allow one to calculate thermodynamic functions and electron transport coefficients 
for such plasmas with an accuracy required for theoretical interpretation of observations. 
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1 Overview 



Envelopes of neutron stars (NSs) are divided in the inner and outer envelopes. Proper- 
ties of the inner envelopes are outlined in a companion paper M. Here we focus on the 
outer ones, at typical densities p ~ 10 4 — 10 11 gem -3 and temperatures T ~ 10 4 — 10 9 K. 



These envelopes are relatively thin, but they affect significantly NS evolution. In par- 



ticular, they provide thermal insulation of stellar interiors. 



Bulk properties of the NS envelopes (pressure P, internal energy U) are determined 
mainly by degenerate electron gas. The electrons are relativistic at x T > 1, where 
x T = hkp/mc ~ (p^Z/A) 1 / 3 , kp = (Sn 2 ^) 1 / 3 is the Fermi wave number, m is the 
electron mass, n e = Zn- X is the electron number density, p 6 = p/10 6 gem -3 , Z and A 
are the ion charge and mass numbers, respectively. Degeneracy is strong at T Tp = 



(eF — mc 2 )/k-Q, where e-p = mc 2 ^ r and 7 r = ^JT+x 2 is the Lorentz factor. 

The ionic component of the plasma can form either strongly coupled Coulomb liquid 
or solid, classical or quantum. The classical ion coupling parameter is Y = [Ze] 2 j ak^T, 
where a is the ion-sphere radius (4^a 3 = ). The ions form a crystal if V exceeds some 
critical value T m (see below) . The quantization of ionic motion becomes important at 
T ^ T p , where T p = h y / 47m i Z 2 'e 2 / 'M /k^ is the ion plasma temperature, and M is 
the ion mass. The ionic contribution determines specific heat Cy, unless T <^T p . 

The ions and electrons are coupled together through the electron response (screen- 
ing). The inverse screening (Thomas-Fermi) length is /ctf — (47re 2 dn e / dp) 1 ! 2 , where \i 
is the electron chemical potential. If T T F , then p m ep and k^p/kp m 4(a/ir) 7 r /x r , 
where a is the fine structure constant. Note that 7 r /a; r — > 1 at x x ^> 1; therefore the 
screening effects do not vanish even at high densities. 

Thermodynamic properties are altered by the magnetic field B in the case where 
the Landau quantization of transverse electron motion is important. The field is called 
strongly quantizing when it sets all electrons on the ground Landau level. This occurs 

at p < p B and T < T B (see, e.g., ref. §), where p B « 7045 (A/Z) B 3 ^ 2 gcm- 3 ,T B w 
1.343 x 10 8 ( J Bi 2 /7r) K, and B 12 = B/l¥ 2 G. Weakly quantizing fields (p B < p) do not 
significantly alter P and U but cause oscillations of Cy, other second-order quantities, 
and electron transport coefficients with increasing density. 
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Characteristic domains in the p-T diagram are shown in Figure |l]. The short- 
dashed lines on the left panel indicate the region of partial ionization. The dot-dashed 
lines correspond to L = 1 (upper lines) and L = 175 (liquid/solid phase transition). 
Long-dashed lines on the right panel separate three p-T regions where the magnetic 
field is strongly quantizing (to the left of ps and considerably below Tb), weakly 
quantizing (to the right of ps at T < Tb), or classical (much above Tb). 

The dotted lines in 
Fig. [I] show profiles T(p) in 
the envelope of a "canoni- 
cal" NS with mass 1.4 M Q 
and radius 10 km, and 
with an effective surface 
temperature T s = 5x 10 5 K 
or 2 x 10 6 K (the values of 
T s6 = T s /10 6 K are marked 
near these curves). Typ- 
ical temperatures of iso- 
lated NSs are believed to 
lie in the hatched region 
between these two lines. In 
a magnetized NS, T(p) de- 
pends on strength as well 
as direction of the mag- 
netic field. Therefore on 
the right panel we show 
two dotted curves for each value of T s6 : the lower curve of each pair corresponds 
to the heat propagation along the field lines (||, - i.e., near the magnetic poles) and 
the upper one to the transverse propagation (_L, near the magnetic equator). 

2 Equation of state 

Thermodynamic functions of the electron gas at arbitrary degeneracy are expressed 
through the well known Fermi-Dirac integrals. For astrophysical use, it is convenient 
to employ analytic fitting formulae for these functions presented, e.g., in ref. j3|. 

Nonideal (exchange and correlation) corrections for nonrelativistic electrons at fi- 
nite temperature have been calculated and parameterized in ref. [Q. For the relativistic 
electrons at low T, an analytic expansion of the exchange corrections is given, e.g., in 
ref. || (in this case, the correlation corrections are negligible). A smooth interpolation 
between these two cases has been constructed in ref. ||. 

For the ionic component at V ^> 1, the main contribution to the thermodynamic 
functions comes from the ion correlations. Strongly coupled one-component Coulomb 
plasmas (OCP) of ions in the uniform electron background have been studied by 
many authors. The thermodynamic functions of the classical OCP liquid have been 
calculated [7], || at T > 1 and parameterized 0, [6[] for < T < 200. The latter 
parameterization ensures accuracy ~ 10~ 3 /cbT (per particle). 

For the classical Coulomb crystal, accurate numerical results and fitting formulae 
to the free energy F (with anharmonic corrections taken into account) have been 
presented in refs. pTCfl . A comparison of the latter results with the fit || for the 
liquid yields the liquid-solid phase transition at L m = 175.0 ±0.4 (curiously, this value 
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Fig. 1: Characteristic p-T domains for outer NS envelopes 
composed of iron. Left panel: non-magnetic plasma; right 
panel: B = 10 12 G. Shown are /9-dependences of Tp, T p , Tb, 
Pb (see text); lines where F = 1 and P = 175 (dot-dashed 
curves); lines where effective Z equals 15 or 20 (short-dashed 
lines on the left panel) . Dotted lines and hatched regions show 
temperature profiles in NSs at different cooling stages and (on 
the right panel) at different inclinations of the magnetic field. 
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could be derived in ref. @, but it was first noticed in ref. f|). 

In the solid phase, quantum effects may be important at the considered temper- 
atures. These effects can be most easily taken into account in the approximation of 
a harmonic Coulomb crystal. A convenient analytic approximation is provided by a 



model [|TTJ which treats two phonon modes as Debye modes and the third one as an 
Einstein mode. Numerical calculations [12] of the harmonic-lattice contributions to 
F, U, Cy, and the mean-square ion displacement are reproduced by this model at 
arbitrary T/T p within several percent, which is sufficient for many applications. 

Finally, there is a contribution from 
ion-electron (ie) correlations, which 
can be treated as polarization of the 
"jellium" of the degenerate electrons. 
In a Coulomb liquid, this contribu- 
tion has 
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been evaluated in refs. 
14] using a perturbation theory which 
involves the electron dielectric func- 
tion and the static structure factor of 
ions. This approximation is justified 
at A>tf *C fcp (which is the case at 
P6 ^ 1CT 2 ). At smaller densities the 
ie contribution in the Coulomb liquid 
has been calculated in ref. using the 
HNC technique (heavy dots in Fig. |2|). 

In a Coulomb solid, the ie contribu- 
tion has been evaluated in re f. || using 
the perturbation theory [|13|, [l4f and a 
model lattice structure factor. These 
results, though approximate, indicate 
that the polarization corrections may 
be rather important. In particular, T m 
may be shifted by ~ 10% when the ie 
corrections in the liquid and solid phases are taken into account. 

The numerical results have been fitted |3|, by analytic functions of x T and T. The 
approximation || for the OCP liquid is shown in Fig. ^] by solid lines. Since the ie 
correction is not dominant, the accuracy of the fitting formulae (< 10%) is sufficient 
for most of astrophysical applications. An alternative Pade approximation, proposed 
in ref. \TE\ for the electron-ion fluid, is also accurate at x r < 1 for small and large T, 
but it is less accurate at intermediate Y and inapplicable at x r > 1 (dashed lines). 

As mentioned above, quantizing magnetic fields significantly affect the equation of 
state. In the fully ionized dense plasma, the main effects are described by the model 
of ideal electron gas, studied by many authors (e.g., |2|, |T6|1). 
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Fig. 2: Nonideal ie contribution to the free en- 
ergy of carbon (normalized to N\(Ze) 2 /a, where 
iV"i is the number of ions) at 8 values of x r marked 
near the curves. HNC results [|| and relativis- 
ts perturbation theory |Q (YS) are compared 
with analytic fits |J| (SB) and @ (PC). 



3 Electron conductivities 

Heat conduction in NS envelopes is provided mainly by the electrons. The thermal 
conductivity tensor k is determined mainly by the electron-ion scattering. Calculation 
of k should take into account specific features of the Coulomb plasmas in NS envelopes, 
quite different from the terrestrial liquid and solid metals: (i) in the liquid phase, 
an incipient long-range order emerges at T > 10 2 ; (ii) in the solid phase, there are 
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usually many Brillouin zones within the Fermi surface; (iii) in the solid phase, the 
approximation of one-phonon scattering fails near the melting, which necessitates the 
use of a more general expression for the structure factor of ions. These features have 



been taken into account in ref. fl7|| , where analytic approximations to the electron 
transport coefficients have also been derived. 

Magnetic field hampers electron transport across the field lines: transverse ther- 
mal and electrical conductivities are reduced by orders of magnitude in typical NS 
envelopes. Strongly quantizing field significantly changes also longitudinal conductivi- 
ties; weakly quantizing field causes de Haas-van Alphen oscillations. These effects have 
been outlined, e.g., in ||. A unified treatment of the electron transport coefficients 
in the domains of classical, weakly quantizing, and strongly quantizing magnetic field 
has been developed in A Fortran code which implements this treatment is avail- 
able at http://www.ioffe.rssi.ru/astro/conduct/. Using this code and solving 
the thermal diffusion equation, we have calculated the T(p) profiles shown in Fig. [I]. 
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